function res = fn_inveta(m, h0)

if h0>2e-3
    res = sinh(h0/2)^2/(sinh(h0/2*(m+1))^2 - sinh(h0/2)^2);
else
    res = 1/(m^2 + 2*m);
end
     
end